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Abstract 

We study classical spin networks with group SU2. In the first part, using gaus- 
sian integrals, we compute their generating series in the case where the edges are 
equipped with holonomies; this generalizes Westbury's formula. In the second part, 
we use an integral formula for the square of the spin network and perform stationary 
phase approximation under some non-degeneracy hypothesis. This gives a precise 
asymptotic behavior when the labels are rescaled by a constant going to infinity. 

1 Introduction 

A classical spin network is a pair (F, c) where T is a trivalent graph equipped 
with a cyclic ordering of the edges around vertices and c is a map from the 
edges to the natural numbers called coloring and satisfying some simple 
conditions. According to Penrose |P71j . one may associate to such a pair a 
rational number (r, c) obtained by contracting a tensor with values in some 
representations of SU2. 

When r is a theta graph (Q) or a tetrahedron {(^), these quantities 
were introduced by Racah and Wigner in 1940-1950 for the study of atomic 
spectra |W71j . Later, Ponzano and Regge used it as a discrete model for 
gravity. Immediately, physicists were interested in the study of the asymp- 
totic behavior of the quantity (P, kc) when k goes to infinity. This behavior 
corresponds to a classical limit of quantum mechanics and is expected to be 
related to euclidean geometric quantities. 

In the nineties physicists used spin networks in the spin foam models 
for quantum gravity [B98J and the study of the asymptotical behavior was 
extended from the 3j (thcta) and 6j (tetrahedron) to more complicated net- 
works as the 9j (complete bipartite 3,3 graph) lOj, 15j (skeleta of the 4 and 
5 simplices), see |BS03l IFL0 2] . Mathemat icians got interested in spin net- 
works through the work of Kirillov, Reshetikhin [KR89j and Kauffman [K94j 
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who introduced and studied their "quantum versions" which are a key in- 
gredient in the construction of quantum invariants of knots and 3-manifolds. 
In this paper, we wiU not deal with these "quantum versions" but restrict 
to the classical case. 

The first rigorous proof for the asymptotical behavior of the 6j-symbols 
were obtained (in the so-called euclidean case) by Roberts [Rj and then re- 
obtained and extended using different techniques |C101 IGV09| IWTOSj . For 
general graphs, Garoufalidis and Van der Veen proved that the generating 
series of the sequence k i— )• (F, kc) is a G-function, implying that the sequence 
(F, kc) is of Nilsson type and thus that the asymptotic behavior does exist 
|GV09j ■ Abdesselam obtained estimations on the growth of spin- network 
evaluations, specially for generalized drum graphs, see |A09j . 

Another interesting approach to the study of spin networks was proposed 
by Westbury who computed in |W98j the generating series of spin networks 
as follows (we refer to Section [2] for notation). Let Rx = C[[Xci,a € A]] be 
the ring of formal variables associated to the angles of F. For any coloring c, 
we denote by Ce, Cq and X*^ respectively the color of the edge e, that of the an- 
gle a and the monomial Haeyi^a"- Then, writing ((F,c)) = (F,c)i=p^^^, 
we define the following formal series: 

z{r) = Y,{{r,c))x^. 

c 

Given (5 C F, a subgraph which is a disjoint union of cycles, we denote by cs 
the coloring which associates to an angle 1 or if 5 contains the two edges 

forming the angle or not respectively. Let then Pr = Yl ■ 

<5cr 

Theorem (Westbury, |W98j ). If T is a planar graph then Z{T) = P^'^ ■ 

This was generalized to all trivalent graphs by Garoufalidis and Van der 
Veen |GV09j . Here below we describe our results which may be presented 
in two independent parts. 



1.1 Generating series with holonomies 

In gauge theory it is natural to consider spin networks whose edges are 
decorated by a holonomy ip in SL2(C): we shall denote it by (F,c, ■0) (see 



Definition 2.3) and its renormalization by {{T,c,ip)) (defined as above). In 
the present paper we generalize Westbury's theorem to the case of a general 
graph F equipped with any holonomy ■0 and we give a closed formula for 
the generating series Z{T,tp) = ^^((F, c, 

Fix a numbering of the vertices of F compatible with the planar structure, 
see Figure [2] and let F' be the graph obtained from F by blowing up vertices, 
i.e. replacing ^vL^"^^^^- Vertices of F' are in 1-to-l correspondence with the 
set H of half-edges of F and the holonomy may be seen as a map ip : H ^ 
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SL2(C). Then, we set = for every half-edge and for any pair of half- 
edges g, h, we define bg^h '■ Fg x ^ Chy the formula hg^h[zg,Wg, z^, Wh) = 

i{ZgWh - ZhWg). 

We define the following i?x-valued quadratic form on F/i : 
Q{x) = 2 hg^rM-^Xg,i,^^XH)X^ + 2 ^ hg 

a:g—^h e:g-^h 

The expression a : g ^ h means that a is an angle between the half-edges g 
and h such that the opposite vertex of g is lower than the opposite vertex of 
h. In the same way, e : g ^ h means that the edge e contains the half-edges 
g and h in such a way that the vertex contained in g is lower than the vertex 
contained in h. Using Gaussian integrals, we prove the following result: 

Theorem 1.1. Given a planar graph T with holonomy ip : H ^ SL2(C), 
the following formula holds Z(r,'i/') = det(Q)~^/^. 

In this expression, the determinant is computed in the canonical ba- 
sis. We remark that strictly speaking, the formula makes sense only for 
holonomies in SL2(M) because of the indeterminacy in the square root. By 
analytic continuation, the formula holds in general. The non-planar case, 
can be easily treated as a corollary of Theorem |l.lt see Subsection 3.3 As a 



corollary we observe that Z(T, V')"'^ is a polynomial with integer coefficients 
in the entries oi ip. In Theorem |3.1[ we provide a combinatorial interpreta- 
tion of our formula which extends Westbury's formula in the case when ^lJ 
takes values in the diagonal subgroup of SL2(C). 

1.2 Asymptotics from integral formulas 

In his book |W71j . Wigner showed that the square of a 6j-symbol may be 
computed by a simple integral formula over 4 copies of SU2 . Barrett observed 
that this formula may be generalized to any graph pBS03j. In Section[4| after 
recalling the integral formula, we compute the generating series of squares 
of spin networks: 

Theorem 1.2. Let (r,c, ^) be a spin network equipped with a holonomy 
with values in SL2(C). Define 

(r,c,'0)^ 



[r,c,v^] 



v:{ei,e2,e3) 



where by v : (£1,62,63) we indicate the edges touching v and (0, Ce^, Cej, Ceg) 
is the value of the Q- graph colored by Ce^ , Cgj , ■ Considering the generating 
series W{V,%I:) = X^ct^^c, V"] V", the following holds in Ry = C[Ye,e G E]: 

dg 



T^(r,v) 



e={hi,h2)={v,w) 
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Using Kirillov trace formula, the integral formulas may be transformed 
in order to apply the stationary phase approximation. This method was 
applied by |BS031 IFL02j to compute the asymptotical behavior of some spin 
networks but it faces some technical difficulties because of the existence of 



so-called "degenerate configurations" . In Section |1.3| we propose a different 
transformation which allows us to treat uniformly all configurations corre- 
sponding to critical points. We describe precisely the critical points of the 
integrand and compute the associated Hessian. Under suitable genericity 



hypotheses on F and c described in Section 5.2, we compute the dominat- 
ing terms in the asymptotical development of [F, kc, 1] for general F. Fix a 
trivalent graph F and a coloring c : — >• N. 

Let I be the set of maps P from oriented edges of F to 5"^ which satisfy 
the following relations: 

- Denoting by — e the edge e with opposite orientation, we have P_e = 

-Pe- 

- For all vertex v with outgoing edges ei, 62, 63 we have Yli CeiPe^ = 0. 
Given P £ I we define rp(C) = Ce\\Pe x for ^ E ^^^^ ^^(^^^) ^ 

T.e:{v,w)(^e\\Pe X {^v - ^wW for (^^) £ ^^^yR^. In thesc formulas, we 
denoted by x the bracket in the cross-product in M^. 

Given a pair (P, Q) of non-isometric elements of /, we define its phase 
iction T in 



function r in Subsection 5.2 as a map from E to S'V{±1}- Then we set for 

p3. 



e ^ 

Then, given a quadratic form q on M", we denote by det'(g) the determinant 
of the restriction of the quadratic form to the orthogonal of the kernel of q, 
that is the product of all non-zero eigenvalues of the matrix of q. We also 
set det'(gp,Q) = limK^i(K - l)-^ det'(g^Q). 

Theorem 1.3. Let (F,c) be a colored graph satisfying the conditions of 



Subsection 5.2. Then denoting by N the opposite of the Euler characteristic 



ofV, one has 

_ (2iV)3/^ /^det(rp)V^ 



(^A;3)A^-iV/^^det'(gp)i/2 



(P,(9)e/2/±i,P^Q 

The above formula was numerically tested in the case of the tetrahedron 
and of course its results are in line with the known asymptotics ([R], |GV09j . 
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|BS03] . |C10j ). The non-degeneracy conditions of Subsection 5.2 are equiv- 
alent to the non-vanishing of all determinants in the formula of Theorem 
Ol 



- The quantity det(rp) is zero only if the configuration P is planar which 
can occur only for very special values of c. 

- The non-vanishing of det'(gp) is equivalent to the infinitesimal rigidity 
of the configuration P. In particular, the hypothesis do not hold if the 
set / is not discrete: this happens for instance for the regular cube 
or more generally, the cube whose edges are colored by the lengths of 
Bricard's flexible octahedron. 

- We do not have a geometric interpretation of the determinant det' (qp^q) 
but in our numerical experiments on the spin-networks formed by the 
1-skeleton of a tetrahedron and of a triangular prism this determinant 
was non-zero. In any case, the conditions define a Zariski open set of 
configurations. 

We believe that describing when the non-degeneracy conditions hold is a 
very difficult task as it contains the problem of the fiexibility of polyhedra, 
a notoriously hard problem. Still, we expect that for planar graphs whose 
colors correspond to the lengths of a generic convex configuration of the 
dual graph, there is a simple geometric condition ensuring that the non- 
degeneracy conditions hold but this question is not addressed in this article. 



1.3 Some questions 

- The generating series Z{T, ip)'^ has coefficients which are integral poly- 
nomials in the entries of ifj. Is it still true for Z{r,ij)) as suggested by 
the abelian case? 

- Is there a combinatorial interpretation of Z{r, -0) for general ipl 

- Find a direct relation between the series Z{T,ip) and W{T,ip)7 

- Find sufficient conditions on F and c which ensure that the non-degeneracy 
hypothesis hold, in particular in the case of convex polyhedra. 

- Interpret geometrically the leading terms in the asymptotic formula of 
Theorem II. 31 
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2 Equivalent descriptions of spin network evalua- 
tions 

Let r be a trivalent graph with possibly loops and multiple edges. We 
denote by E the set of edges, V the set of vertices and we divide each edge 
e € E into two subarcs called half-edges. We will then let H be the set of 
half-edges commonly described as pairs (e, v) where e is an edge and v is 
an end of e. We shall moreover orient a priori each half-edge h = (e, v) so 
that it goes out of the vertex v. We assume that for each vertex the set of 
half-edges incoming to that vertex have a cyclic order. We call angle a pair 
of half-edges touching the same vertex and denote by A the set of angles. 
In all the article, we will suppose that F is connected and contains at least 
a vertex; we will denote by N the opposite of the Euler characteristic of T: 
we have #V = 2N, #E = 3N and #H = #A = 6N. 

We will denote by admissible coloring a map c : E ^ N satisfying the 
triangle conditions: 

Vu : {i, j, k), Ci + Cj -I- Cfc G 2N and Cj < Cj + (T) 

For convenience, we wrote v : {i,j,k) meaning that the edges i,j,k are in- 
coming at V with that cyclic order. Associated to each coloring c there is 
an internal coloring, i.e. a map from A to N also denoted by c and defined 
as follows: if a is the angle between edges i and j around v then we set 

_ ci+c^ cfc ^ remark than one can recover the original coloring from 
the internal coloring: to avoid confusion, we will denote edges with latin 
letters and angles with greek ones. 

Let iV,uj) be a complex symplectic vector space of rank 2 and SL(y) be 
its symmetry group. A discrete connection is a map i/; : H SL{V); we 
define the gauge group as the group of maps {g : V JI E ^ SL(T^)}. An 
element g of the gauge group acts on a discrete connection -0 as follows: 
for each h = (e,v) & H, (g ■ = getphQv^- Two discrete connections 
01,-02 '■ E SL(y) are said to be gauge equivalent if they are in the same 
orbit of the gauge group. 

Definition 2.1 (Holonomy). - A holonomy on F is an equivalence class 
of discrete connections on F. 

- The trivial holonomy (denoted by 1) is the class of the constant discrete 
connection defined by -0^ = 1 G SL(y) for all half-edges h G H. 

- Let 7 be an oriented path in F described as a sequence of half-edges 
hi, ... , hn- We define the holonomy -0(7) of along 7 to be the product 
V"^" • • • V'^^ where £j is 1 or -1 if the half-edge is oriented coherently with 
7 or not. If V'l and 02 are gauge equivalent and 7 is closed then -01 (7) 
and V'2(7) are conjugate. 
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Figure 1: Invariant tensor at vertices with a = 3,6 = 5,c = 4. 

In the fohowing section, we give two equivalent descriptions of the spin 
network evaluation which consists in associating to a triple (F, c, ■0) a com- 
plex number (F, c, tp). We will refer to these descriptions as the abstract and 
computational descriptions. 

2.1 Abstract description 

Let (y, u) be a complex symplectic vector space of rank 2. We will consider 
that V is an odd superspace. We refer to Appendix for a basic review of 
supersymmetry adapted to our purposes. 

The symplectic form is considered as a supersymmetric map ui : V0V ^ 
C in the sense that u o 99(12) = For any integer n £ N, uj induces a map 
w®" : V^"" ^ C defined by the formula 

n 

Uj'^"'{vi (g • • • (g) -Unjtfl <X' • • • (8) Wn) = JJ „,+l_i , Ifj) . 

1=1 

Denote by Vn the subspace of V^" consisting in anti-supersymmetric 
tensors (or, equivalently, symmetric in the standard sense); the form w®" 
restricts to a supersymmetric form Wn : y„ (g 14 — )• C. The vector space 
Vn (of parity n) is the (n + l)-dimensional irreducible representation of the 
group SL(y): we will sometimes use the notation pn ■ SL(V) — )• End{Vn)- 
Denote by (jJ~^ the unique element of V g) V such that the contraction of 
the two middle terms in a; g> is the identity of V. 

Let a,b,c be three integers. It is well known that the set of SL(y)- 
invariant elements in 14 g) 14 g) is 1-dimensional if a, b, c satisfy the triangle 
conditions (T) unless it is 0. One can find an explicit generator £a,b,c S 

g) Vb g" 14 given by the symmetrization of the element 

UJa,b,c = (a;-l)®(''+^-^)/2 ^ (^-l)®(6+c-a)/2 ^ (^-l)®{a+c-6)/2 

where the supersymmetric tensor product is reordered as in Figure [l] Re- 
mark that the sign of the permutation which reorders the a + b + c factors 
in this tensor product is 

Then we define £a,b,c — (Lla ® lib ® ^c)'^a,b,c where IIq is the anti- 
supersymmetrization (i.e. standard symmetrization) map projecting V®"- 
onto Va- 



7 



One may check that this element is supersymmetric as if wc permute 
cychcally a, b and c to the left the result is multiplied by (—1)" = (— 
which is the sign of the cycle 14 (8) H ® K ^ H ® K ® according to the 
supersymmetric rule. 

Definition 2.2 (Spin Network). Let F be a trivalent graph and c : E 

be an admissible coloring. Then, (F, c, 1) is the result of the supersymmetric 

contraction of 

((^Wcg, £ci,Cj,Ck)- 

e&E veV,v:{i,j,k) 

By supersymmetric contraction, we mean that we reorder the tensors on 
the right hand side according to the sign rule so that factors corresponding 
to the same edge are consecutive, and then contract with the maps a;„. 

Given a discrete connection tp : H Sh{V), let us define the following 
endomorphism Hol^ by the formula: 

(H) Pci {'^i,v)^Pcj {'<Pj,v)^Pck e End (g) Vc^ ® K,- <8 Kfc 

veV,v:{i,j,k) \veV,v:{i,j,k) J 

Definition 2.3 (Spin Network with holonomy). Let F be a trivalent graph, 
tp : H ^ SL(y) be a discrete connection and c :£'—)• N be an admissi- 
ble coloring. Then, (F, c, tp) is the result of the following supersymmetric 
contraction: 

((g)a;ce,Holv,( (g) ec„c,,c,)). 

eSE veV,v:{i,j,k) 

One can check directly that this definition does not depend on the gauge 
equivalence class of ijj. 



2.2 Computational description 

The computational description follows directly from the abstract one, taking 
care of the signs. 

Set y = and uj{v, w) = det{v, w). Then V* = and we denote by z 
and w the corresponding coordinates. In this way, uj~^ eVi^iV corresponds 
to the linear polynomial on V* x V* which in coordinates {zi,wi, Z2,W2) 
reads as Z2W1 — ziW2- 

The form a;„ : (g) V„ — C is a linear form on the space of homogeneous 
polynomials of bidegree (n, n) on V* xV*. Here we identify homogeneous 
polynomials and symmetric tensors via: 

n! ^-^ 
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Figure 2: Planar presentation of the tetrahedron 



where X^, = z \i k < i and w otherwise). One checks directly that in 
coordinates, u;„ is expressed by: 

\ d d d d 

^ __(_____ \n 

" n!^ dzi dw2 dz2 dwi 

In the same way, the element ea,fe,c S Va Vb (d) Vc corresponds to a 
polynomial on V* x V* x V* which in coordinates (zi, wi, 2:2, 1^)2, 23, W3) 
reads as 

£a,b,c = {Z2WI - Z1W2) ^2 (Z3W2 - Z2W3) ^2 (Z3W1 - Z1W3) "^2 . 

Finally, let us compute the sign needed to perform the final contraction. 
Suppose that the trivalent graph is presented as in Figure [2] this presenta- 
tion induces an ordering of the set of half-edges H. For any edge e G E, let 
ei and 62 be the respective left and right half-edges of e and for any ver- 
tex vinV, let vi,V2,V3 be the three half-edges incoming to v in increasing 
order. Writing in line the tensor product of invariant elements ea,6,c give 
a big tensor which has to be reordered in such a way that half edges are 
packed in pairs. Let C be the set of pairs of crossing edges: the sign of this 
permutation is (-l)S{e,<='}ex ^ec,, _ 

We compute then: 

(r,c,i) = (-i)^i.^'}.x'^^'^^'Y\—(—— ^^V' 

Ce\^ \dZei dWe2 dZe^dWe^J 
{z^^Wy^ - Zv^W^^) ^ {Zv3Wv2 - Zy^Wy^) 2 (1) 



{Z-D^W-D-^ Z^-^^W-u^) 2 

In order to use the machinery of Gaussian integration, it will be more 
comfortable to introduce i = \f—\ in our formulas. The formula is un- 
changed if we replace Wn and ea,b,c by tJ„ and ea,6,c where we set 

_ i"- d d d d 
n!^ dzi dw2 dz2 dwi 

a-|-b-j-c a^h — c b~t-c — a a^c — h 

£a,b,c = i 2 (^ziW2 - Z2W1) 2 (^Z2W3 - Z3W2) 2 [ziW^ - Z3W1) 2 . 
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The formula for (F, c, ifj) is obtained from Equation ([T| by replacing the 
variables {zi^Wi) at the vertex v with 'ip~^{zi,Wi). 

3 Generating series via gaussian integrals 



m 



In Subsection |3.1| we will use gaussian integration to prove Theorem 1.1 
the case of a planar graph F. This is a generalization of Westbury's result to 
the case when graphs are equipped with holonomies. In Subsections |3.2| we 
will provide a topological interpretation of the theorem when the holonomy 
is diagonal and recover Westbury's original result (always in the case of 
planar graphs). In Subsection 3.3 we will extend Theorem [lT] to the case of 
non-planar graphs. 

3.1 Computing the generating series of a spin network 

Let r be a planar trivalent graph that we present as in Figure [2 It may be 
easily seen that the sign correction due to crossings in Equation (ij) is 1. Let 
Ffi be a copy of associated to each half-edge h, we will always consider 
the standard density on these spaces and then omit it in the notation, and 
let F = 0^ F/j. For two half-edges g, h we set bg^h ■ Fg x Fh ^ C as being 
given in coordinates by the expression i{zgWh — z^Wg) and h~\ : F* xF^ ^ C 
as —i(-d_J^ 8 d \ 

^ dzg dwfi dzfi dwg ' ' 

Define quadratic forms P on F* and Q on F respectively by P = 
—2 Yli ^gh ^^"^ <5 = 2 ^ Xabg^h- The notation a : g ^ h and e : g ^ h 

e:g—>h a:g—^h 

means that a and e are composed of the half-edges g, h which appear in 
that order from left to right in Figure [2] For the moment, Xa should be 
interpreted as a real parameter. Note that the quadratic form on F is 
expressed by = 2 6g,h- 

Consider P as a differential operator P°p on C°°(F, C) and develop 
(exp(ip)°^'exp(iQ))|o (see Appendix [b|) . Collecting monomials in the vari- 
ables Xa, by Equation ([T]) one sees that the coefficient of fla ^oT is ((F, c)) 
(defined in the Introduction). Therefore we have: 

Z(F,l) = (exp(^P)°^'exp(^g))|o. 

Consider a deformation of Q which is non-degenerate and has positive 
real part. For concreteness, we can pick Qq = ^^'^ ~ Q''^ 

eQo. We define = ( exp(^P)°P exp(|Q^)) |o so that Z(F, 1) = lime^o^e- 
Then replacing Q by and ^P hy exp(2P), Formula ([l3]) gives 

Z, = (27r)-"/2det(g-i)^/2^^ exp(^P(x) - ^Q-Hx))dx. 
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We apply now Formula (12) to the integral, noting that the quadratic 
form — P is still non-degenerate and has positive real part. Hence, we 
have 

' detiQe^-P) ^ ' 

Letting e go to 0, we find that Z{T, 1) = det(Qo - Q-P)"^/^. 

Suppose that ip is represented by a discrete connection on T with values 
in SL2(M). By Formula ([T]), we know how to adapt the construction: for each 
angle a connecting two half-edges g and h, we need to replace hg^h '■ FgXFh ^ 
C by bg^h{\(j~^Xg,ipJ^^Xh). We denote by the resulting quadratic form. 
By the assumption that ip lives in SL2(M), P takes again only imaginary 
values. Hence, the argument above repeats exactly and we obtain Z(T, ip) = 
dei{Qo-Q^P)-^/'^. The general case, that is, for ip taking values in SL2(C) 
follows by analytic continuation. 

One can simplify this formula by remarking that the matrix of P in the 
canonical basis satisfy P~^ = —P and moreover, we have det(P) = 1. We 
obtain the formula of Theorem 11.11 



Theorem ( 1.1 ). 

Z(r,V) = det(P + Q^)-i/2. 

3.2 A generalization of Westbury's theorem to the case of 
diagonal holonomies 

Suppose that the holonomy ^p can be represented by a connection with 
values in the subgroup DcSL2(M) of diagonal matrices, we introduce a map 

i : i7 — )■ M* such that for all h, one has 'Ph = i^!! ) in the basis {zh,Wh). 

In this case, we can extend Westbury's theorem [T] as follows. 

Let C(r) be the set of all oriented curves immersed in F which pass 
over an edge of F either 0, 1 or 2 times, in the latter case with opposite 
orientations. Given 7 G C(F) we denote by cr(7) the number of crossings 
modulo 2 of the corresponding immersion. 

Let X"' = Haci (hence each angle may appear 0, 1 or 2 times) and 

1^(7) = n/iC7 th^"''^^ i^ili ^) being 1 if 7 crosses h in the positive direction, 
— 1 otherwise). The following result generalizes Westbury's result to the case 
of holonomies with values in diagonal matrices: 

Theorem 3.1. Let T be a planar graph equipped with an abelian holonomy 
ip. Then: 

Z(F,V) = ( (-l)"(^)Tr(7)X^(c))-' 
7ec(r) 



Theorem |3.1| descends directly from Corollary C.2 and Proposition C.6 
of Appendix [C] where we also show how to recover Westbury's theorem from 
it. 
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3.3 The non-planar case 

Let now F be a non-planar graph and let us fix a diagram of T as in Figure 
[2] containing crossings xi,...Xk- By the supersymmetric rules, for each 
coloring c of F, a crossing Xj between edges ei and 62 induces a factor 
^_X)CeiCe2 wjiicii we denote by s{xi,c). The previous results allow us to 
compute the "wrong" generating series for F where the signs coming from 
crossings are not taken into account, that is we can compute: WiTjip) = 
^^((F, c, V'))-'^'^ s(xj,c). To fix these signs, use the identity 

(-1)'''' = ^(1 + (-1)'^ + (-1)' - i-ir+'), Va, b G Z. 

More explicitly, for each edge e let ae,/3e be the leftmost angles formed by 
e and let Opg : C[[X]] — t- C[[X]] be the automorphism that changes the 
signs of Xa^ and Xjs^ . Then, for a crossing x between edges ei and 62 define 
: C[[X]] ^ C[[X]] as 5^ = ^(Id+Op,^ + Op,^ - Op,^ o Op J. Then we 
recover Z{T,^p) = '^c{{r,c,ip))X'^ as 

Z(F,V) = 5.i0...o5,,(T^(F,^)). 

In particular, when ■(/; is trivial this recovers Garoufalidis and Van der Veen's 
extension to non-planar graphs of Westbury's theorem |GV09] . 

4 An integral formula for the square of spin net- 
work evaluations and its analysis 

In this section we deduce an integral formula for the square of the spin- 
network evaluation. In the case of the tetrahedron, this formula was already 
known to Wigner (see |W71| ) . It was generalized by Barrett in [BS03j to any 
graph. For questions of clarity and normalizations, we provide a proof of it, 
then we derive a formula for the corresponding generating series. From now 
on we will no longer work in the super-symmetric category: being interested 
in squares of spin- network evaluations, sign matters are irrelevant. 

4.1 Derivation of the integral formula 

Let (y,u;,h) be a symplectic complex and hermitian vector space of rank 
2. Both structures are supposed to be compatible in the sense that there 
is a hermitian basis (61,62) of V such that a;(ei,e2) = 1. We also suppose 
that the hermitian product is anti-linear on the left and denote by SU(y) 
the symmetry group of the whole structure. Naturally, y®"- has a hermitian 
structure as its subspace Vn ■ Denote also by h the hermitian product induced 
on the spaces Vn and their tensor products. 

We would like to replace contractions with the forms cOn by scalar prod- 
ucts with some vectors. The compatibility of the hermitian structure and 
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the symplectic form give immediately the formula: h{uj w) = uj{w, v) 

or h{w,v w) = (jj{v,w) where we set vo = —u:~^. Let Wn = (n„ 
nn)(g)I Li^i,2n+i-i) where n„ : V®"^ — )■ T^®"' is the projector on symmet- 
ric tensors in the non-supersymmetric sense. After some computation, we 
obtain for any v,w ^ Vn, h{zun,v w) = Univ, w). 

Then we can compute the spin network evaluation using scalar product 
instead of contraction, more precisely we have up to sign: 



{T,c) = h\(^wc^, ec„c„cj. (2) 

eS-B veV,v:(i,j,k) 

Denote by G the group SU(y) and for any integer n, we denote by pn '■ 
G — >• U(y„) the induced representation of G. We also write dg for the Haar 
measure on G satisfying J^dg = 1. 
Consider the element 

'P= Pn{g)® Pn{g)dg £^nd{Vn®Vn). 

Jg 

We check directly that V satisfies V o V = V and V* = V. Moreover, 
the image of V is the one dimensional space {Vn ® Ki)*^- Hence, V is the 
orthogonal projector on Ctij^. Using th.6 formulci hi^zUf^^ 
we have h{wn-, Wn) = {'^n-, ^n)- But in the contraction of ujn ® one can 
first contract the middle terms. One find n times the evaluation of uj w 
which is equal to — Idy*. Then, the middle contraction is (— l)"Idy^. The 
final contraction computes the super-trace. Hence, one find h{wn-,'^n) = 
n+1. Finally, we obtained that for any v G Vn^Vn, 'Piv) = u)ci7". 

In the same spirit, consider for any triple (a, b, c) satisfying the triangular 
relations the element 

S= pa{g)(E> pb{g)(^ pc{g)dg eEnd{Va(E>Vb(E>Vc). 
Jg 

The same reasoning as above shows that for any v G Va <Si Vb <^ Vc one has 

(a,6,c) 

Lemma 4.1. 



£{v) = where we set (a,6, c) = h{ea,b,c-,£a,b,c)- 



(a,6,c) = |(e,(a,6,c))| 



^ a+b+c _|_ a+h—c | ^ a—h+c -j | ^ —a+b+c ^ j 



Proof. A proof of it may be found for instance in |BL81j or computed from 
the generating series given by Theorem |W98j for the theta graph, taking 
care of the renormalization from ((0, (a, h, c))) to (0, (a, 6, c)). □ 

In equation Q, replace e hy £ and w hy V and consider the scalar 
product h {^^^^Ve:®.u^y £v\- usiug the fact that all matrices involved 
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are self-adjoint, one sees that it equals Tr(0ggg7'e o ^^gy^i))- For any 
edge e colored by Cg, let (e) = Cg + 1 and Ve be the projector acting on 
Vc^ ® Vc^. Similarly, for any vertex v whose incoming edges arc i.j. k, we set 
(v) = {ci,Cj,Ck) and denote by £y the projector acting on ® Vc^ (g) Vc^. 

Let [r,c] =ne(e)'IV((g),ggPeO(8),g^f.) =ne(e)^((8)eei?^e,(g)„gyf.). 

The first term in the scalar product defining [F, c] is the orthogonal 
projector on ^^w^" and the second term is the orthogonal projector on 

£v hence, we have: 

|r,^l=n>r- (3) 

Here we use the fact that h induces a hermitian form on End(W) = WCg) W* 
and that the scalar product of projectors Pu and on unit vectors u and 
V is |(u, On the other hand, we can invert hermitian product and 

integration finding 



v:(i,j,k) 



[r,c] = n(e)/ l[ Tr:^,ig^ge)dg= [ H TT,Sgv9;n^)dg. 

Remark that the last equation holds even if c is a non admissible coloring 
on r (in which case both sides of the equations are 0) . In order to prove the 
last equality, we used the following lemma: 

Lemma 4.2. For any integer c and all elements A,BeG one has 

f TrdAg) 'Itc{Bg)dg = TrdAP-^). 
JG C-h i 

Proof. We have 

/ Tic{Ag)T:v^{Bg)dg = f Trc{AB-^g)Trc{g)dg 
JG JG 

= Trpc (^AB-^ p^{g) Trc{g)dg'^ . 

Writing U = pc{g) Tii:c{g)dg one obtains that U commutes with G hence 
U is proportional to the identity. Moreover, Tr{U) = jQTrc{g)^dg = 1 
hence U = ^ Idy^ which proves the lemma. □ 

As a summary we obtained the following formula, noting that (F, c) is 
real: 

(r,c)2 = n(^) n Tr,A9v9^')d9 = ll{^)[r,c] 
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Notice that we need to orient arbitrarily the edges of F to write this formula: 
this detail will be important in the next section. 

Let R = C[[Ye,e G E]] be the ring of formal variables associated to the 
edges. To any coloring c, we denote by the monomial Hee-B ■ For any 
A in G, the following identity holds: 

Consider then VF(r) = X^cl-'^' '']^^ where the sum is taken also over 
non-admissible colorings by setting [T, c] = in that case. We deduce the 
following formula: 



W(T) = y = / (4) 

cad^i^sibie n.(^) n det(id-yw) 

e:{v,w) 



Example 4.3. Consider the case of the graph Q. Then for any coloring 
c given by three colors a, b, c, one has by construction (©, a, b, c) = (a, 6, c) 
and hence the generating series Q is simply 



(1 - YiY2){l - l2l3)(l - Y^Y^) Jg nL det(Id -Y,g) ' 



c, admissible 

On the other hand, the coefficient of Yf'Y2Y^ in the above integral is 
Tra(5) Trf,(g) Trc((7)d5(. Orthonormality of the characters in L'^{G,dg) 
and the Clebsch-Gordan rules imply that this integral is 1 if a, b, c are ad- 
missible and if not. 

We can generalize the formula to the holonomy case: fix a holonomy 
tjj : H ^ G. Then, one has to replace at each vertex v the tensor eci,cj,ck 
with Hol^(eci,c,-,cJ = Pci{ipi,v) <^ Pc,{ipj,v) <^ Pcki'4'k,v)£ci,cj,Ck- Then in the 
computation above, we need to replace the projector £ with the projector 
on Hol^ (eci,cj,Cfc) that is Hol^^'Hol^"^. Therefore we have: 

[T,c,^p] = T[{e) / rr Trc,{geTphgviJh^)dg 

e n={v,e) 



n Trc,(V'e,t;5i>V'e,^^e,»«fi't«Ve,i)d5' 



Since as before the above equalities still hold if c is not admissible (all terms 
are 0), defining l^(r, ^) as in Equation Q, the following holds: 

Theorem 4.4. Let (F, -;/;) he a trivalent graph equipped with a connection 
with values in SL2{C), and let Ry = C[[Ye,e £ E]]. Then the following 
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equality holds in Ry: 



w(r,^) 



dg 

e:{v,w) 



4.2 Comparison between generating series 

Let r be a trivalent graph. We denote by Rx be the sub-algebra of Rx 
hnearly generated by X'^ for admissible colorings c and by Ry the sub- 
algebra of Ry linearly generated by for admissible colorings c. We define 
a morphism (anti-linear on the left) (•, •) : Rx x Rx — >• Ry in the following 
way: if c and c' are distinct then {X^, X^ ) = 0. If not, we set 

yc 

{X^^X") = 



n ((©1 Cej^ , , ££3)) 
i':(ei,e2,e3) 



The equation ^ relating (F, c) and [F, c] generalizes directly to the case 
with holonomies. Summing over all colorings we find: 



(Z(F,^),Z(F,V^)) = Ty(F,V'). (5) 

It would be interesting to interpret this formula at the level of generating 
series for which we should interpret the product (•, •) as an hermitian product 
over Rx given by some integral formula. We do not address this question 
here. 

We end this section with some remarks on the orthogonality of spin- 
network as functions of the holonomy. Let M be the moduli space of 
holonomies on F with values in G = SU2- It is equal to the quotient of 
qH ^qVUE. Haar measure on produces a measure ^ on A4. 

The map from M to C sending ip to (F, c, ip) is a polynomial function 
and it is well-known that the family (F,c, •) form an orthogonal basis of 
L^(A4,/i). The equation (|3]) relating (F,c) and [F,c] generalizes directly to 
the case with holonomies. We can use it to compute the normalization of 
the basis: 



We have jQTTc{AgBg^^)dg = ^^Trc{A)Tic{B) by an argument similar 
to Lemma |4.2[ If we apply this formula to each term of the product, the 
term He^^) S^^^ inverted. By the orthonormality of characters and Clebsh- 
Gordan rule, we find |(F, c, V')pd|u(V') = llv{'^)Ile{(^)~^ ■ 
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5 Asymptotics of spin-network evaluations 



In this section, we compute the first order of the asymptotic behavior of 
[r, kc, 1] (the square of the spin-network evaluation divided by nD(^))' 



general F. In Subsection 5.1 we transform the integral formula to adapt it to 



stationary phase approximation and describe the critical set of the integrand. 



Then we discuss some sufficient non-degeneracy hypotheses (Subsection 5.2). 



The computation of the Hessian occupies Subsection |5.3[ in Subsection 5.4 



we apply the stationary phase method whereas Subsection 5.4.3 ends the 
proof. 

5.1 Description of the critical set 

From now on, to simplify the notation, we will write for any v,w G V : 
h{v,w) = {v,w). Denote by the unit sphere of V and dv the Haar 
measure on it (i.e. such that J^a Idv = 1). 

Lemma 5.1. For any g £ G, Tr„(g) = (n -|- 1) Jg3{v, p{g)v)'^dv. 

Proof. Recall that 0„ Vn is the algebra of polynomials on V* and that G 
acts on it by algebra morphisms. For any v £ , we have {v^,w^) = 
{v,w)"'. Consider the endomorphism U of Vn given hy U = J^s v'^{v'^, •)dv. 
This endomorphism clearly commutes with the action of G, hence it is 
proportional to Idv„- As Tr(t>"(t>"', •)) = 1, we find Tr(?7) = 1 hence 
U = :^ldv„- Now, given any g G G, we have Tr„(g) = Tr pn{g) = 
in + l)TTpn{g)U = {n + l) Jg,{v'', p{g)v'')dv = {n + 1) f^^iv, p{g)v)''dv. □ 

We deduce the following formula: 

[r,c] = TT(e) / / TT {gvUe,gwUey''dgdu. 

The notation e : v ^ w means that e is an oriented edge joining v to w, Ue 
is an element of corresponding to the edge e. 

In order to analyze this integral, we restrict the domain of the integral 
such that the integrand does not vanish and such that symmetries behave 
nicely: we will see that this restriction does not affect the asymptotic anal- 
ysis. More precisely, let X he the following subset of G^ x (5'^)^ and 
TT : 5^ — 7- S"^ be the Hopf fibration. The sphere is either identified to 
the projective space ¥(y) or to the unit sphere of the Lie algebra ^ of G 
(equipped with the scalar product = — ^ Tr(^^) ). 

X = {{gy,Ue) G G^ X (S^)^ such that \/e : v ^ w, {g^Ue, g^Ue) ^ 

and the family (7r(ne))ee£; has rank at least 2 in Q} 
Let F : X — )• C/2i7rZ be the map F{g,u) = Yl Celn{g^Ue, gwUe) and 

write Pe = T^{ue) E ^ if the orientation of e coincides with the chosen one 
and Pe = —7r{ue) if not. We have the following proposition: 
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Proposition 5.2. The critical points of F are the elements {g,u) G X 
satisfying 

1. For all edges e : v ^ w, gvUe = T^gwUe for some £ . 

2. For all vertices v with edges 61,62,63 going out of v, one has Ce^Pei + 

F^2 ~1~ Fq^ — . 

Proof. Derivate F with respect to the variable Ue in the direction Ve- We 
have DF{ve) = i^g^^^^g^ue) iidvVe, gwUe) + {gvUe,gwVe))- Denote by the 
element of satisfying {ue,u-^) = and det(ue,ii^) = 1. Then, write 
g~^gw in the basis {ue,uj:). Taking Ve = uj: and Ve = iujr , we find that 
9v^9w is diagonal in the basis (uejW^)- This is equivalent to the first item 
of the proposition. 

Consider an element ^ in Lie(G) and replace g^ by (71,6*^ for some t > 0. 
Then compute the derivative of F at t = 0. Three terms contribute to 
the sum, corresponding to the edges incoming to v. Suppose that they are 
oriented from v to fi,f2,f3 by edges denoted by 61,62,63 respectively. For 
simplicity, we note g = g^, gi = gvi, Ui = and a = Ce,. 

We have F'(0) = EiiC*^g^ = ELi c^(Cni, n,). The last equality 
comes from the first item. Let Hi be the element of End(y) acting by i on 
Ui and —i on n^. Then, from Tr(^nj) = 2i{ui,^Ui) we have that X^i '^j^i is 
0. Using the standard identification between anti-hermitian operators and 
Q which sends Ilj to Pj, we obtain the second item of the proposition. □ 

Suppose that the coloring c satisfy the strict triangle inequalities (T): 
then at any vertex, the triple Pei , -fe2 j -Pea has rank 2. In particular, the 
corresponding pair {gv,Ue) is in X. We will make this assumption in the 
sequel. 

The integral presents some obvious symmetries: the integrand depends 
on Ue through its image tt{U(,) E S'^. Hence, we can replace the integral 
over (S^)^ by an integral over (S^)^ integrating over the fiber of the Hopf 
fibration (S^)^. Moreover replacing {gv,Ue) by {ggv,Ue) do not change the 
integral as replacing {g^^Ue) by {gvg~^,gue)- Hence, the integral can be 
performed over Y = X/{S^)^ x G x G where {ae,g,h) acts on {gy,Ue) by 
{ggyh~^, h{aeUe))- Notice that the stabilizer of the action of (S^)^ x G x G 
on X is {±1}, hence the quotient y is a smooth manifold of dimension 
12N — 6. Let us denote by F : y — )• C/2i7rZ the induced map. 

We denote by / C /G the isometry classes of tuples (Pe)ees satisfying 



the equation (2) of Proposition 5.2 They encode the critical points of F in 
the following way: 

Proposition 5.3. The set C of critical points of F are in bijection with 
the set of triples {P, Q, (gv)) in I x I x ^G^/{ibl}^ where for all half-edges 

{e,v) we have gvPe = Qe- Moreover, the map C ^ I x I which forgets the 
third term is surjective and its fibers have cardinality 2^^"-^. 
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We remark that if > 1 , I is never reduced to 1 element as if P belongs 
to I, —P is an element of / distinct from P. 

Proof. Let {g, u) be a critical point of F. Then the family Pg = 7r(ue) belongs 
to /. For any half-edge (e,u), set Qe = 7r{gvUe) = 7r{gwUe)- For any vertex 
V with edges 61,62,63 incoming to it, the relation X^jCe^Pei = implies 
X^jCgjQei = 0, hence family (Qe) also belongs to I. Conversely, given two 
families (Pe), (Qe) in /, we choose vectors Ue,Se € satisfying 7r{ue) = Pe 
and 7r(se) = Qe ■ Then for every vertex v surrounded by ei, 62, 63, we need 
to find an element g^ € G such that for all i, g^Ug. is proportional to Sg^- As 
we required that the coefficients Ce^ satisfy the strict triangular inequalities, 
the pair P^^ , P^^ is linearly independent in ^. As Qei , Qe2 have same lengths 
and same angle by hypothesis, there is a unique rotation of G which sends 
Pe; to Qe- for all i. This rotation lifts in two possible ways in G giving 
our choices for g^j. Note that the action of{±l}GGxG divides by 2 the 
number of critical points. □ 



5.2 Non-degeneracy hypothesis 

We would like to understand which conditions ensure that the critical points 

of F are isolated in Y. At first, we concentrate on the family of vectors (Pe) 

associated to critical points. 

Let £ be the set of maps from E to G. We consider it as an Euclidian 

space where the scalar product is induced from the scalar product = 
-5Tr(^2) g T^is 

space will then be identified alternatively as the chain 
or cochain complex of F with coefficients in ^. As edges of E are oriented, 
there is a boundary operator d : £ ^ Co(T , Q) and we set Ti = kerd. We see 
that any critical point in X produces an element C of by setting (^e = CgPe- 
Conversely, any element C, oiT-L with |^e| = Ce produces an element of I (see 



Proposition 5.3). 



A tangent vector to Q is an element Q' ^ % such that (Ce,Ce) = fo'^ 
all 6 (because |Ce| = Ce)- We interpret this formula as the scalar product of 
{Ce^e) and C'- On the other hand, a tangent vector correspond to a global 
symmetry if and only if there exists ^ & Q such that Ce = C ^ Ce for all e. We 
count that the dimension of is 3(A^+1), the number of conditions imposed 
by C is 'iN (one for each edge), whereas the symmetry gives 3 dimensions 
(because the family Ce has rank at least 2). Hence, an element of /G 
associated to P is isolated if the family {Pe^e) projected orthogonally on 
% is linearly independent. With this interpretation in mind, we make the 
following assumption: 

Hypothesis 1. Let (ce) he a coloring ofT such that all triangle inequalities 
in (T) are strict. We suppose that all elements Q ^ % such that |Ce| = Cg 
satisfy the non- degeneracy condition that the family (n-^(Ce^e))eGS is lin- 
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early independent where 11-^ : £ ^ % is the orthogonal projector on % . This 
implies that I is a finite set. 

Remark 5.4. The preceding hypothesis is equivalent to saying that for any 
{^v) G Co(r, Q), if Pe X [iw — ^v) = ^ for aU edges e : v ^ w then all the 
are equal. Let us show it: 

The image of d : Co{T,Q) — )• £ is orthogonal to the space 71. Let U 
be the subspace U = (Be^Pe C S. The non-degeneracy assumption is that 
the orthogonal projection oiU onTi is injective. Hence, by standard linear 
algebra, the orthogonal projection of U-^ on Ti-^ is surjective. This implies 
that if some vector (^g in 'K-'^ satisfies Pg x Ce = for all e then = 0. 

For any critical point x, one can define the phase function : E ^ 
in the following way. Set x = {g^, Ug). Then, by the first item of Proposition 
|5.2[ g^Ue and gwUe are proportional unit vectors, moreover the phase factor 
Te = {dvUe, 9wUe) depends only on the class of the critical point x in Y. We 
remark that the generalized phase of a critical point indexed by (P, Q, gv) S 
I X I X depends on g^ only up to a sign and is necessarily equal to ±1 if 
P = Q. We assume that this is the only case when the phase function takes 
the value ±1: 

Hypothesis 2. For any P,Q, distinct elements of I, the phase function of 
the associated critical points does not take the value ±1. 

Let P, Q be distinct elements of /, and r be the associated phase function. 
Write Te = e*^'=. Then, 6^ is well-defined modulo vr. Consider the following 
quadratic form on : 

q = ^2Ce {-iC0t{9e)\\Qe X {^v - ^w)\\'^ + i{Qe, X ^^)) . 
e 

Hypothesis 3. For any P,Q distinct elements of I, the quadratic form q 
has co-rank 6. 

Let us develop an example which is the main motivation of this section: 

Example 5.5. Suppose that T is a planar graph and that / = {P, — P}. 
Then there is up to isometry a unique polyhedron A C G whose 1-skeleton 
is dual to r (hence whose faces are triangles) and such that for any oriented 
edge e of r the dual edge in A is vectorially equal to CePe- in particular, it 
has length Cg. 

The non-degeneracy condition (HI) is equivalent to an infinitesimal 
rigidity condition on A. For instance, if A is convex, this condition is au- 
tomatically satisfied by Cauchy's theorem (see for instance |AZ09| . Chapter 
13). The generalized phase function has the following nice interpretation: 
let {gv,Ue) be a critical point associated to the pair (P, — P). This means 
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that TT{ue) = Pe and TT^g^Ue) = —Pe- Let w be a vertex of F, dual to a face 
of A. Suppose that v has three outgoing edges 61,62,63. Then, hfts the 
unique rotation mapping P^^ to —Pei, that is the rotation of vr in the plane 
supporting F^. Now, given an edge 6, one has Te = {gyUe, QwUe) = ±e'^ 
where 26 is the angle of the rotation g~^gw around Pe- Hence, 6 is the 
angle modulo vr between the faces F^ and F^. In particular, if A is a non- 
degenerate polyhedron, its corresponding graph satisfies (H2). We do not 
know which geometric condition on A would imply that (H3) is valid. 



5.3 Second variation formula 

Given a point (i^^, Ue) G X, we can build a a coordinate system around it us- 
ing the variables G , (Ae) G and Oe G {S^)^ . The parametrization 
is given by {e^^g^, ae^e(Ae)) where for any u ^ and A G C we set 

= + ^^^) = (1 - 1^1 V2)n + A^.^ + o(|A|2). 

Let = gvUe and notice that gyUe{Xe) = s%{\e)- Let ix^ = {s%.,^yS%), 
H% = {s%-^ ,^yS%) , Te = In what follows, the values oi gy,Ue, are 

to be considered as fixed and ^i,,Ae (and consequently those of 
variable. Notice that G M, that Te was already defined in Proposition [5[2| 
and that F : X ^ R does not depend on the coordinates {oe}, so we set 
them to 1 in the following: 

Proposition 5.6. The Taylor expansion of F in a neighborhood of a point 
{gv,Ue) £ X is up to the second order: 



F(6^"5^,,Ue(Ae)) = ^CeluTe ^Cei(a 

e e 

+ Y,Ce-(i4 + Xe){fil + Xe) 

- Y,Ce{\\e\^ + \tff/2 + \^fJ^/2 + Ye^xl + KK) 

e 

V e 

In particular, on critical points J^e'^'i^v ~ ^w) = ^'^^ ^ = '^e'^- 

e ■ e 

Ml) ~1Xy 

moreover since (gf,,g f„) = Te, then ((s^)"*", (•s^)"'") = ( = Tg^^ on critical 
points by Proposition 5.2 ). To compute the first terms of F{e^^ g^, Ue{\e)) = 
^ Cg ln(e^"s^(Ae), 6^™s^(Ae)), we first compute 
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Then we get (e^"s^(Ae 



(Ae))=re(l- 



I A, 



-ie.iV2- 



/x^Ae) + Te{n'i, + Ae)(Mw + ^e) + /i-o.t. Taking the logarithm, 



expanding it up to order 2 terms, and recalling that 
(x^)^ + we get the formula of the proposition. The last statement is a 
consequence of Proposition |5.2[ □ 



To obtain a nicer formula, in the case of a critical point, we introduce the 
variables zf, = ^% and = T~'^fi'^. After some computation, by Proposition 



5.6, one has that on critical points Iiess{F)x{(,, X) = —q{S,,X) where: 



Ke,A) = ^Ce(2(l-r-2)|Aep + 2Aer-2(zS,-z,-) + 2Ae(z:-z^) 



+ \Zv 



1 -r; 



(6) 



+- 



Te + r: 



z;r + z^zS,-4^^) if^e /IVeG^. (7) 



Tf. - Tr 



-1 1^10 



5.4 Applying the stationary phase method 

Let us now perform the stationary phase approximation when replacing the 
coloring c by kc and letting k go to infinity in Jy exp(kF)dfi. In this formula, 
the measure ^ is obtained from the Haar measure on (S^)^^ x (S^)^^ by 
integration over the action of (S^)^ x G x G/{ibl}, equipped with its Haar 
density. In this subsection as well as in the following we rely on the notation 
and machinery recalled in Appendix [Bj 

By assumption, any critical point x is isolated in Y. Provided that 
Hess{F)x is non-degenerate, we can apply the stationary phase expansion 
theorem (Theorem 7.7.5 of [H83]) to F which is a smooth function with 
non-positive real part, and obtain that the local contribution of x to [T, kc] 
when k goes to infinity is: 



The proof of Theorem |1.3| consists then in computing for any critical 
point the gaussian integral I{q) where q is the opposite of the Hessian com- 
puted in the preceding subsection. We observe that in Equation ^ this 
quadratic form is expressed in terms of {X^, z^,z^} and that its restriction 
to {z^ = z^ = 0} is diagonal. Thus we would like to perform a partial 
integration on the Ae terms which form local coordinates for the (5^)^-part 

of y. 

There are two different cases to handle: in the first case one has Tg = 1 
for all e in The quadratic form q is then degenerate with respect to the 
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A coordinates and we deform it with a parameter k. In the second case, one 
has Tg 7^ 1 for ah e, thanks to Hypothesis 2. Fix a critical point described 



by a pair {P,Q) G I x I (as in Proposition 5.3); to avoid a cumbersome 



notation, in what follows we will suppress the indices P, Q except in the 



definition of the quadratic forms used in the statement of Theorem 1.3 



5.4.1 Deformation and degeneracies 

The kernel of q contains the image of ^ x ^. Let and r] be the two 
variables acting on {gv,Ue) by {e^gve~^ ,e^Uu)- For an edge e : v ^ w set 
Pe = {u-^,r]Ue) and qe = (s^"'-,^s^). The infinitesimal action of ^ and r] in 
the tangent space at {gv,Ue) is — gvVdv^^v)- O^^s computes in terms of 
Pe and Qe this action as Ae = Pe, = = (s^-^, - g^rig~'^)s%) = qe - Pe 
and z^ = qe- r-'^Pe- 

Consider the deformation q^ of q obtained by replacing Tg by = uTe- 
One recovers the original one by letting n G (l,+oo) go to 1. After some 
computations one may check that we have qK\gxg = 4(k — ^)Yle'^e\Pe\'^ + 
o{k — 1) Hence, denoting by rp the quadratic form on Q given by rp{rf) = 

Cg|pgp we can remove the indeterminacy in rj as explained in Appendix 
E We have: 



^1 1{qK\gxg) I{rp) 
5.4.2 Partial integration 

We now remark that this deformation also allows to integrate partially 
over the coordinates (Ag). More precisely, let 0"^ be the restriction of q^ 
to the A-coordinates, setting cJk = X^g 2cg(l — r^^)|Agp, one obtains = 
I{aK)I{q'^) where 

o'^-Vc /^^^^^^^^tl^Ue _ |2 , e^_-e 

y .Z—/ \ —l\^V ^w\ ' ^V^W ^V^ 

e \TeK, TeK 

We can reformulate the preceding quadratic form by introducing only 
operations on Q. For instance, one has = \Qe x ^^p. We also have 
^v^w ~ ^w^v = ~^ T^(J^Cv(,w) = '^i{Qe, X where H is the matrix whose 
eigenvalues in the basis {sf,,sf,-^) are i and —i and x S,w = [£,v,£,w]- This 
may be proved by direct computation in the basis (5^,5^-*-). Hence, we can 
write: 



1 



1P,Q = ( y^^"-i \Qe X {iv - + 2i(Qg,e, Xi^)]. (9) 
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First case: = 1 for all e in E. 
Write = ^^^Q + q' where we set q' = 2ice{Qe, x ^w) and 

|2 



Using remark 5.4 the kernel of qq is the subspace Q given by the 
equations = i for all v. When k goes to 1, the quadratic form q' be- 
comes negligible. Hence, we have ~ (k - lf^-^/'^I{qQ), and usmg 
= /(|Ap,/ig^^)/7r = 2 one computes: 

(4k -4)3/2 Ijqg) _ 2^'-'^I{qQ) 

I(rp) (4k - 4)3^ neCe(K- 1)3/2-3^ ^(^p) lie ' 

Second case: 7^ 1 for all e in 

Setting n = 1 in q'^, we get the quadratic form q": 

ip,Q = (-^cot(ee)|nQx(c. - e»)P + 2i(ge,e. X u) ■ (10) 



e 



By hypothesis (H3) this form has corank 6 hence, the following formula 
makes sense: 



3Ar 



(4k -4)3/2 (4k -4)3/2 2 
/ = hm ^ z — i I(ai)I(q'^) = hm ^ — ^ ^/(g 



^.^ 23(K-l)3/2/(g'^) _ ^.^ 23-3^i^ne^e(/^-l)^/'J(g") 
«™1 /(rp) He Ce(l - Te"^) ~ Ue Sm{9e)I{rp) 



5.4.3 Collecting the critical points 

In this section, we apply the previous computations to all critical points 
and collect the results. There are two cases: whether the critical point 
corresponds to a pair (P, P) or to a pair (P, Q) with P ^ Q. 

First case: Fix a SA'^-tuple {ue) representing P^.. All critical points asso- 
ciated to the pair (P, P) have the form ((—1)^", Ue) where G {if}- In the 
preceding subsection, we obtained that the contribution of a single critical 
point (1, Ue) to ([s]) is 2^~^^ I (qp) / 1 (rp) , remarking that the term Yie can- 
cels and using the fact that F{1, Ue) = 0. The other pairs ((— 1)^", Ue) differ 
only by the value of e^^ = W {^—\){<^v+ew)kce _ Adding the contributions 



and dividing by 2 (see Proposition |5.3[ ), we get: 

11 ^ ^ ; Urp) Urp) ' 

£i> i;:(ei, 62,63) 

We note that the contribution of (— P, —P) is the same because q-p = 
qp and r„p = rp. We will multiply the result by 2 when summing over 

P e //{±i}. 
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Second case: In the case when P ^ Q, we obtained in Equation (11) 
that the contribution of a critical point {gv,Ue) is 

23-3iV^^(^_l)3/2j( . ) 

K^i I{rp)ll^sm{ee) 

As before, taking into account all critical points associated to the same pair 
{P,Q) amounts in multiplying the result by 2^^~^. 

When replacing {P,Q) by {—P,—Q), we also have r_p = rp. Given 
a critical point {gy^Ue) corresponding to {P,Q), one checks directly that 
the point {gy,u^) corresponds to {—P,—Q) and that the phase function 
gets inverted. Letting qpQ be defined as in Equation one sees that 
q'^P^Q = —QpQ and since qpQ is purely imaginary when k — )• 1, it holds: 
lini(K — l)^/^/(— (7pq) = lini(K — l)^/^/(gpg). Collecting the contributions 
of {P, Q) and (— -P, — Q), we get the following formula: 

23-^Re(i^ lim(K - l)3/2j(g« ^)eiE.(fcce+i)e, 



/(rp)nesm(0e) 

It remains to compute explicitely the Gaussian integrals involved in the 
preceding computation. Considering on Q x Q the usual Haar measure, one 
has to divide the result by 2 because of the isotropy subgroup {±1} £ GxG. 
Using Appendix |B| we can replace all occurrences of I by a usual det in the 
euclidean basis of G and its tensor powers. Denoting by det' the product of 
the non-zero eigenvalues of a matrix, we have 

I{rp) = (2/7r)i/2^et(rp)-^/2 

liqp) = (2iV)3/2(2/7r)^-i/2det'(gp)-V2 

/(gp,Q,M) = (2iV)3/2(2A)^-i/2det'((?^,Q)-i/2 

where the factors {2N)'^^'^ are due to the fact that the Haar density on the 
kernels of qp and of qpQ is fJ-Haar = (^2N)^/'^ because each vector of the kernel 
is represented in Q"^^ by 2A^-copies of the same vector of Q. 



A Super-symmetric rules 

We refer to |DM96| for a detailed discussion. We collect here some definitions 
and warnings important in the article. A super vector space F is a direct sum 
of two finite dimensional complex vector spaces Vq and Vi called respectively 
even and odd part. We will always suppose that our spaces are homogeneous 
that is one of the two components vanishes, in particular any element x € 
V has a degree |x| which is or 1 depending on the parity of V. Super 
vector spaces form a category where morphism are linear maps respecting 
the decomposition. 
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The tensor produet V W is defined as the vector space whose even 
part is Vb ® Wo e ¥l (gi Wi and odd part is Vi (^Wq ®Vo 0Wi. There is an 
isomorphism cv,w ■ V ^ W (E>V sending x y to (— (g) x. The 
even vector space C is neutral for tensor product. 

This allows to define the unordered tensor product of homogeneous super 
spaces in the following way: Given a finite family of vector spaces {Vi)i^i of 
parity pi, we set 

(g) = lun ((^(^,)-ioa : Va ^ K')- 

The projective system is defined as follows: let n be the cardinality of 
/, for any bijection a : {1, . . . , n} — > /, set Vo- = 0^=1 Then, define an 
isomorphism V'((t')-1o<t • Kr' by the formula 

frivi (8) • • • (8) fn) = (— 1)*'"ti <Si ■■ ■<SiVr„ where s = ^ PiPii. 

i<i' ,Ti>T^t 

There is an internal functor Horn satisfying the following adjunction 
formula: 

Hom(^, Hom fy, W)) = Hom(?7 ® V, W) 

The space Hom(i7, V) is the space of all linear maps where a map is consid- 
ered even if it respects the parity and odd if it reverses it. In particular, if 
we set U* = Hom([/, C) we have an isomorphism Hom([/, V) = V U* and 
an evaluation map cvu : U* (Si U C 

With these identifications, bilinear forms are elements of Hom(C/ 0V, C) = 
Hom(C/, V*) = V* i^U*. The inversion of the terms should be noticed. 

There is another tricky issue: the natural isomorphism 9u ■ U ^ (U*)* 
is not the identity but sends x to (—1)1^1 With that convention, we have the 
identity evu* =ev[/ o cu,u* ° (^[/^ ® 1). Unformally, this formula is explained 
by the transposition of the terms in the Gelfand transform: 9{x){X) = 
(-l)l^ll^lA(x). 

B Gaussian integrals and densities 

In this subsection let F be a real vector space of dimension n, F* be its 
dual, Q a quadratic form on F and S{F) the symmetric algebra of F. If 
Q is non-degenerate, one can identify F and F* and thus transport it to a 
quadratic form denoted on F*. Note that if F is equipped with a basis 
and F* with the dual, the matrices expressing Q and in these bases are 
inverse to each other. 

Definition B.l (Densities). Let B{F) be the set of bases of F. A density 
/X on F is a map /x : B{F) M such that for each A G GL{F) it holds 
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fi{A • 6) = I det{A)\^{b). The set of densities is a real 1-dimensional vector 
space denoted |A|(F). A density on a manifold M is a continuous section of 
the bundle \A\{TM). 

Fixing b G B{F), defines an isomorphism F ~ and thus F can be 
equipped with a density which we denote iJ,euc which satisfies fieucib) = 1. 
If F is equipped with a density fi and Q : F ^ C is a quadratic form, 
then one defines det(Q,/i) G C as det((5(5))//i(6)^ where Q{b) is the matrix 
expressing Q in the basis b oi F (this clearly does not depend on b). 

Lemma B.2 (Gaussian Integrals). Given a density fi on F and a quadratic 
form Q such that Re[Q) > 0, then 

I{Q) = [ exp(-^)d/x = Z^""^"^' (12) 

where the square root is the analytical extension of the positive one on the 
set of real and positive quadratic forms. 

An element P of S{V*) may be interpreted either as a polynomial func- 
tion on y or as a differential operator on C^{V*;C). To distinguish the 
two cases we shall denote P°p the element P interpreted as a differential 
operator. The following is an ubiquitous generalization of the Gaussian 
integration formula (IT2|): 



Proposition B.3 (Fourier transforms of Gaussian functions). 

/ P(x)exp(--Q(x))dM= \ J (poPexp{-Q-'))\o (13) 

If K C F is a, A;-dimensional subspace one has |A(F)| = |A(K)| (g)]^ 
\A{F/K)\: indeed given bases bx and bp/^^ for K and F/K one can construct 
easily a basis of 6 of F using the inclusion of K and choosing an arbitrary 
complement to K (the choice of the complement affects b only up to elements 
of SL{F)). Thus given densities fi and /xj? on F and K respectively, the 
quotient density Hf/k on F/K is defined so that /i = fix ® fip/K- 

Suppose that the quadratic form Q is degenerate. Then denoting by K 



its kernel, we can apply the formula (12) to the reduced quadratic form Q 



on F/K. To do that, we need to fix a density on K. Then, we can set 

IiQ,I^K) = fp/j^e^p{-lQ{x))dfiF/K{x)- 

This new integral can be computed without considering the quotient by 
the following perturbative argument. Let Q' be a quadratic form on F with 
positive real part and such that it is non-degenerate on K . Then, for e 
positive and small enough, the quadratic form Q + eQ' is non-degenerate on 
K and on F. Moreover we have 
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I{Q + eQ') 



As a consequence , if a € M is positive, it holds I{aQ,^k) = Oi^2^ I {Q , Hk) ■ 
Suppose now that F = Fi(BF2, that Q is non-degenerate on Fi and that both 
F and Fi are equipped with densities fi and /ii (consequently F2 inherits a 
density /i2)- Let Qi be the restriction of Q to Fi and let A : F2 ^ Fi he 
defined by Q{x,y) = Qi{x, Ay),\/y S F2,yx G Fi. Then, setting Q'{y) = 
Q{y,y)-QiAy,Ay) it holds: 

HQ) = i{Qi)i{Q') (14) 

Given a manifold M equipped with a density /i one defines Jj^j fdfi G 
C, V/ G C°°(M) in the natural way. In particular, a compact Lie group G 
can be equipped by the G-invariant density fiH defined by /^^ IdfiH = 1. If 
a G acts freely on a manifold M equipped with a density ^ preserved by G, 
one defines naturally a density fJ-M/G the quotient. 

Example B.4 (The euclidean density on ^). On Q equipped with the scalar 
product = -i Tr(^2)^ 

f^euc be the density whose value on a orthonor- 
mal basis is 1. Then fi^uc coincides with the density induced by the iden- 
tification of Q with the tangent space to S"^ C at a point and the Haar 
density is fin = 

Example B.5 (The Hopf fibration). Let =SU(2) be the unit sphere in 

and let be the quotient by the diagonal action of S^. Let nf^ 

and IjL^ be the Haar densities on SU(2) and the quotient density on 

respectively. Let also ^euc'^^ , /^euc be the euclidean densities. Then /x^^^^ = 

1 , su(2) , pi _ 1 pi 
2:;^H'euc anu — ^fJ-euc- 



C Around the proof of Theorem 3.1 



The following are general well known facts which we shall apply to interpret 
topologically some of the determinants we will be dealing with. 



C.l Combinatorial interpretation of determinants 

Let G°^ be a graph whose vertices vi, . . .Vn are connected by oriented edges 
dj : Vi — ?■ Vj whose weights are the entries rriij of a n x n matrix M (the 
diagonal terms correspond then to loops). The following is a standard well- 
known fact: 

Lemma C.l. det(M) = (— 1)" X]c(~l)'^'^^('') where c runs over all the 
oriented curves embedded in G°''' and touching each vertex exactly once, w{c) 
is the product of the weights ruij of the oriented edges in c and #c is the 
number of connected components of c. 
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1p2 














Figure 3: Blowing up with holonomies 

Let now G be the graph whose vertices vi, . . .Vn are connected by ex- 
actly one unoriented edge eij and suppose that M is a matrix such that 
ma = Vi. Given a connected oriented curve c embedded in G connecting 



vertices Vi-^ 



let wic) 



' m 



and if c is a 



disconnected union of disjoint oriented curves ci, . . . Ck, let w{c) = Yl^ w{ci). 
Also let a dimer be disjoint union of edges d = ei 
w{d) 



■ m 



U • • • U Ci^j^, and let 
and #d = k. By a configuration of curves 
and dimers we will from now on mean a disjoint union c U d of dimers and 
oriented curves embedded in G such that each vertex is contained exactly in 
one component of c U d; its weight will be w{c U d) = w{c)w{d); similarly a 
configuration of dimers will be a configuration of curves and dimers contain- 
ing no curves. We shall denote Conf(G) the set of configurations of curves 
and dimers on G and DConf(G) the set of dimer configurations. Then the 
following holds: 



Corollary C.2. det(M) 



■l)"EcU.eConf(G)(-l)*^+*'^(c)^(d) 



Finally if M* = -M then det(M) = Pfaff (M)^ and Pfaff (M) can be 
interpreted as counting the dimer configurations in G (see |K63j ): 



Theorem C.3. Pf{M) = J2 



deDCon 



f{G) ±V\'^id)\ 



In the above theorem, the square root is due only to our definition of 
w{d) while the choice of the signs is in general a delicate matter (see |K63j ): 
in our specific cases it will be quite easy to determine it. 



C.2 Proof of Theorem l3Jl 

Let r be a planar spin network equipped with a holonomy with values in the 
diagonal matrices of SL2C and let us use the notation introduced in Subsec- 



tion 3.2 Theorem 3.1 descends directly from Corollary C.2 and Proposition 
IC.6t here below the details. 

Let r' be the graph obtained from F by blowing up the vertices; each 
curve 7 C F can be lifted in a natural way to one in F' (which we will 
keep calling 7) and the connection ^ may be lifted to -0 on F' so that the 
holonomy on a curve and its lift coincide (see Figure [3]) . 
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For any half-edge h of F, the space has a standard basis whose corre- 
sponding coordinates are Zh,Wh- We wih say that a basis element has type 
z OT w. Denote by W the matrix of j{P + Q^) in this basis. Remark that, 
since ip is diagonal, the coefficient Wij does not vanish only if i,j corre- 
spond to adjacent vertices of F' with distinct types. Recall that the set of 
half-edges is ordered and put first the type z basis, then the type w basis. 
Then W is then antidiagonal by blocks. Denote by W^,W'^ the matrices 
indexed by H defined respectively by = Wzg^w^ and W^/j = Wwg,Zh- 
We have det(Ty) = det{W^) det{W^). 

Using the presentation of F' in induced by that of F (Figure [2]) we 
see that the edges corresponding to adjacent angles in F form an angle of 
degrees (so a cusp) pointing in the upwards direction in F' (see Figure [s]); 
we will call such edges internal and the other ones external. 

Remark C.4. If e is an edge of F', let le,re be the half-edges respectively 
at the left and right endpoint of e. If e is an external edge of F' oriented 
from left to right, then the entry of Wi corresponding to e is 1 and that 
corresponding to e equipped with the opposite orientation is —1: indeed both 
entries come from the matrix expressing —iP. Similarly, if a is an internal 
edge of F' oriented from left to right, then the entry of Wi corresponding 
to a is {t'i^^tr^)Xa. and that corresponding to the opposite orientation is 
— {ti^t~J^)Xa: both entries come from the matrix expressing —iQ^. 

Let us denote th = tf^^ . This extends to an automorphism of the ring 
of coeffic ients. We check that W'^ = -W^ = (W^f hence, det{W'^) = 
det(VF^) = det(iy^) because there is an even number of half-edges, so, by 
Theorem to prove Theorem 3.1 it is sufficient to interpret det(l^^) in 



terms of traces of curves. To this purpose, remark that W has on the 



diagonal and we can apply Corollary C.2 with G = F' and M = W^. 



Now remark that if a vertex v of F' is contained in a curve c embedded 
in F' then either c contains an internal and an external edge containing v or 
it contains two internal edges containing v : we will call the latter kind of 
vertices "cusps" of c; also a dimer covering an internal edge will be called 
internal. 

Proposition C.5. Let cL) d be a configuration of curves and dimers con- 
taining each vertex of F' exactly once. Then: 

1. There exists a unique set L of disjoint arcs in F' \ c, whose boundary 
vertices coincide with the set of cusps of c. 

2. Let d' = d\L: the set of internal dimers of d' can be joined by external 
edges of F' to form a unique embedded (possibly disconnected) curve 
c'cF'\(cUL). 

Proof. The first statement is proved by remarking that if e is an edge of 
F' whose endpoint is a cusp of a configuration, then its other endpoint can 
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only be either another cusp (in which case set li = e) or must be contained 
in an internal dimer. By iterating this argument for the edge at the other 
endpoint of that internal dimer one eventually constructs a path li which 
must end at another cusp of c. The second statement is proved by a similar 
inspection. □ 
There is a natural map vr : C(T) — t- Conf(r'): define 7r(7) as cU d where 
c is the oriented curve in T' formed by the edges and angles of T contained 
in 7 exactly once and d is the dimer formed by all the angles contained twice 
in 7 and all the edges not contained in 7. 

Proposition C.6. The map n is a bijection and letting 7r(7) = cD d, it 
holds: 

{-l)#^+*'^w{cUd) = Tr(7)X^ 
where cr{'y) is the number of crossings of j. 

Proof. We construct the reciprocal map s : Conf(r') — )• C{T) by associating 
to a configuration cUd the unique oriented curve 7 which passes once exactly 
over the edges contained in c (and with the same orientation as c), twice 



over all the edges contained in L (see Proposition C.5 applied to cU d) and 
twice over all the edges contained in c' (in that case we consider the curve 
without crossings). This proves our first claim. 

Set it{'j) = cL) d then for any connected component 6 of 7, we define 
w{5) = Y\k=o ^hfe /ife+i where 5 visits the half-edges Hq, . . . ,hn, h^+i = Hq in 
that order. We first claim that the following formula holds: 

{-l)#^+*dw{c)w{d) = {-l)*^w{j). 



We check this formula by looking first at the terms, then the th terms 
and then the signs. 

- The monomial in Xq, occuring in the right hand side is X^, where each 
angle appears as many times it is visited by 7. On the left hand side, 
each time an internal edge a is visited by c produces a monomial Xa 
whereas each time a dimer occupies an internal edge Xa produces a 
monomial X^- As internal edges occupied by a dimer are visited twice 
by 7, these monomials coincide. 



By Remark C.4, if an internal edge is visited twice in opposite direc- 
tions, the resulting monomial in th is equal to 1. Hence, monomials 
occur only for internal edges visited once. One sees that the set of 
internal edges visited once by 7 and cL) d coincide (with orientation) 
hence the resulting monomial is the same. 



By Remark C.4, for any individual dimer do one has w{do) = — 1. 
Hence, one has (— l)*'^w(d) = 1 (we drop here the X^ and t^ parts 
which have already been taken care of). In order to show {—l)'^'^'w{'y) = 
(— l)#'^ti;(c), we recall that to build 7 from c, one needs to add the 
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double arcs L and the double curve c'. For the latter, one sees that the 
weight is 1 because the number of curves is even and each curve goes 
through a even number of edges. Now, any added double arc contains 
exactly one more external edge than internal ones and thus produces a 
-1 factor. On the other hand, either it merges two components of c or 
it splits one component in two hence the formula is proven. 

Let us show finally that {-l)*^w{-i) = {-If'^^'^ Tt{-i)X^ . We only have 
to check the sign, and one can suppose that 7 is connected. Then one has 
to show that Sign w{-i) = {-lf<"i')+^ = (_l)™i°d(7) ^^ere wind(7) is the 
winding number of 7 (here we use the equality wind(7) = cr(7) + 1 mod 2). 
The lift of 7 to an oriented curve immersed in V meets alternatively an 



internal edge and an external one: by Remark C.4 if both are covered in the 
same direction (e.g right to left) then the contribution to ^(7) has positive 
coefficient and if not, then the contribution is negative. On the other hand, 
the tangent vector to 7 rotates either of or ib27r depending on the same 
condition. Hence, one has Sign w{'^) = (— and thus the formula of 
the proposition is proved. □ 

C.3 Recovering Westbury's theorem 



Theorem [T] can be recovered from Theorem 3.1 by remarking that if 



ly he H then the equalities VF^ = -W^ and = {W^ imply {W^f 
-W^ and so det(l^^) = FfaS{W^f. Hence we must show that Pfaff (W^^) 
Y^ccrllacc^a- By Theorem [03] it holds 



pfaff(Tyi)= Yl ^Vwm= Yl ±n^" 

deDConf(r') (ieDConf(r') aCd 

and it is straightforward to check that for each d G DConf(r') the mono- 
mial ^/\w{d)\ equals ±w{c{d)) for a unique (possibly empty or disconnected) 
curve c{d) embedded in F; similarly V7 C T it exists exactly one d G 
DConf (r') such that c{d) = 7. So we have Pfaff (VF^) = J2^^t ^(t) Uac^t 
(with 5(7) = ±1) and we are left to show that 5(7) = 1 V7 C T. Let us first 



remark that by Theorem 3.1 the coefficient of ^(71 U • • • U7fc) = naC7 '^'^ 
det(VF-'^) is 2^. We now use the equafity det(W^^) = Pfaff(M^^)^ and argue 
by induction on A:, li k = 1 the coefficient of w{'y) in det(M^^) is 2 = 25(7) 
and so s{-y) = 1. Now suppose 5(7) = 1 for all 7 C F such that #7 < k; 
then given 7 = 71 U • • • U 7fc+i, the coefficient of ^(7) in det(W^"^) is: 

2^=+! = 2.(7) + 2 Y ^(U ^^)^(U = 2(^7) + I {2''-' - 2)) 

IC{l,...k+l} iel i^I 

and so 5(7) = 1. 



32 



References 



[A09] A. Abdesselam, On the volume conjecture for classical spin networks. 
|arXiv:0904.1734l (2009). 

[AZ09] M. Aigner, G. Ziegler, Proofs from THE BOOK (4tli ed.). Berlin, 
New York: Springer- Verlag (2009). 

[B98] J. Baez, Spin foam models. Class. Quantum Grav. 15 (1998) 1827- 
1858. 

[BS03] J. W. Barrett, C. M. Steele, Asymptotics of relativistic spin net- 
works. Class. Quantum Grav. 20 (2003) 1341-1361. 

[BL81] L. C. Biedenharn, J. D. Louck, Angular momentum in quantum 
physics. Vol. 8 of Encyclopedia of Math. Sci., Cambridge Univ. Press, 
1981. 

[CIO] L. Charles, On the quantization of polygon spaces. Asian J. Math. 14 
(2010), no. 1, 109-152. 

[DM96] P. Deligne, J. W. Morgan, Notes on supersymmetry (following J. 
Bernstein). Quantum fields and strings. A course for mathematicians. 
Vol. 1, 2 (Princeton, NJ, 1996/1997), 41-97, Amer. Math. Soc, Provi- 
dence, RI, 1999. 

[K63] P. W. Kasteleyn, Dimer statistics and phase transitions. Journal of 
Math. Phys. 4, 287-293, No. 2 (1963). 

[FL02] L. Freidel, D. Louapre, Asymptotics of 6j and lOj symbols. 
Class.Quant.Grav. 20 (2003) 1267-1294. 

[GV09] S. Garoufalidis, R. Van der Veen, Asymptotics of classical spin net- 
works, arXiv:0902.3113| (2009), 28 pages. 

[H83] L. Hormander, The Analysis of Linear Partial Differential Operators, 
I. A series of comprehensive Studies in Mathematics, 256, Springer- 
Verlag 1983. 

[K94] L. H. Kauffman, S. L. Lins, Temperley-Lieb recoupling theory and in- 
variants of 3-manifolds. Annals of Mathematics Studies, 134. Princeton 
University Press, Princeton, NJ, 1994. 

[KR89] N. A. Kirillov, N. Y. Reshetikhin, Representations of the algebra 
Uq(sl(2)),q-orthogonal polynomials and invariants of links. Infinite- 
dimensional Lie algebras and groups (Luminy-Marseille, 1988), 285-339, 
Adv. Ser. Math. Phys., 7, World Sci. Publ., Teaneck, NJ, 1989. 

[P71] R. Penrose, Angular momentum: an approach to combinatorial space- 
time, in Quantum Theory and Beyond, ed. T. Bastin, Cambridge Uni- 
versity Press, Cambridge, 1971. 

[R] J. Roberts, Classical 6j-symbols and the tetrahedron. Geom. Topol. 3 
(1999), 21-66. 



33 



[W98] B. Westbury, A generating function for spin network evaluations, 
Knot Theory (Warsaw, 1995), 447-456, Banach Center Publications, 
42, Polish Acad. Sci, Warsaw, 1998. 

[W71] E. P. Wigner, Group Theory and its application to the quantum 
mechanics of atomic spectra. New- York, Academic Press, 1971. 

[WT05] Y. U. Taylor, C. T. Woodward, 6j symbols for Uq(sl2) and non- 
Euclidean tetrahedra. Selecta Math. (N.S.) 11 (2005), no. 3-4, 539-571. 



34 



